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We extract the isovector tensor nucleon form factors, which play an important role in understand- 
ing the transverse spin structure of the nucleon when related to the quark helicity-flip generalized 
parton distributions via their first moments. We employ the light-cone QCD sum rules to leading 
order in QCD and include distribution amplitudes up to twist 6 in order to calculate the three tensor 
form factors Ht, Et and Ht- Our results agree well with those from other approaches in the low 
and high momentum-transfer regions. 
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At the twist- two level the quark structure of the nucleon can be described by three quark distribution functions ff{x), 

O ■ g1{x) and h1[x), where q is the corresponding quark. The unpolarized distribution, ffix), and the spin-dependent 

^^ , helicity distribution, gfix), give a measure of quark contribution to longitudinal momentum and the net helicity of a 

nucleon. These two chiral-even distribution functions are well-known quantities, which can be extracted from inclusive 

deep- inelastic scattering data. The transversity distribution, h\ (x) , measures the net density of quarks having parallel 

polarization to that of the transversely polarized nucleon. h\ {x) , being a chiral-odd quantity, is difficult to measure 

i-^ ' experimentally because it requires a chiral-odd probe to couple in the cross section. The transversely polarized 

Y^. Drell-Yan process [l|, l2| and electro- and photo-production of mesons off the polarized nucleons Q are suggested 

O i' as appropriate ways to measure h'l{x). Recently, h\{x) has been extracted for the first time by an analysis of the 

O , experimental data from Belle [J], HERMES Q and COMPASS [y] collaborations on azimuthal single spin asymmetry 

;^, ' in semi-inclusive deep- inelastic scattering Q • 

At the leading twist, a more general description of the nucleon structure is given by eight generalized parton distribu- 
PsJ ■ tions (GPD): Two chiral-even spin- independent distributions, H and E, two chiral-even spin-dependent distributions, 
^ ] Hi and E^ and four chiral-odd spin-dependent distributions, H!^, E!^, H^ and £'|, Q. Important information about 
the quark-gluon structure of the nucleon is encoded in these quantities. They describe, e.g., how partons arc dis- 
tributed in the transverse plane with respect to motion of the nucleon or the contribution of quark orbital angular 
momentum to total spin of the nucleon (see, e.g., the reviews in Refs. Q)- 

Taking the moments of the GPDs one can obtain the nucleon form factors of the local currents. The first moments of 
the helicity-conserving GPDs, H^ , E^ and H^ , E^ , are constrained by the nucleon form factors of the electromagnetic 
and axial currents, respectively. Similarly, the helicity-flip GPDs, H^, E'^, H^ and E'^, can be related to the nucleon 
form factors of the tensor current by taking the first moment. More explicitly, we have 

/ dxX{x,^,q^)=X{q^), 

C^ I for X = H^, E^ and H^. Note that the first moment of E^(x,(^,q^) vanishes therefore we have E'^(q^) = jlOj . 
Here q = p' — p is the momentum transfer and ^ = —n ■ q/2 denotes the longitudinal momentum transfer, where n is 
a light-like vector. The nucleon matrix elements of the tensor isovector current are parameterized in terms of above 
form factors as follows: 
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where T^i^ = uia ij,,jU — dia ^j^yd is the isovector tensor current, P = p' +p, m^ is the nucleon mass, and X^ ^ = X'"" — X"^ 
for any of the form factors, X = Ht, Et or Ht (hereafter, the / = 1 superscript will be dropped since we will be 
dealing with the I — \ case exclusively). The tensor form factors at Q^ = — g^ = give what are called the tensor 
charges. Previously, they have been investigated in the framework of QCD sum rules [111. |l2|, relativistic quark 
models [ij], axial- vector dominance model [ij] and chiral-quark soliton model(xQSM) [l^. The dependence of the 
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tensor form factors on Q^ have been extracted up to 1 GeV^ from xQSM [lg| and up to a few GeV^ from lattice 
QCD [Hill. 

Our aim in this work is to study the nucleon tensor form factors in the framework of hght-cone QCD sum rules 
method (LCSR) [l^ . LCSR have proved to be rather successful in extracting the values of the hadron form factors 
at high-momentum transfers, e.i?., the electromagnetic and the axial form factors of the nucleon have been calculated 
to leading order and with higher-twist corrections in Refs. 120,l2l|. To our knowledge, while various estimations of 
Ht{Q^) at low- and high- momentum transfers exist, with this study it is the first time in the literature Et{Q^) and 
Ht{Q^) are extracted. 

In this method one starts with the following two-point correlation function: 



<,(p,g) =» / d4xe'«-(0|T[ryAr(0)r^.(x)]|7V(p)), 



(2) 

where rj^^x) is the nucleon interpolating field. There are several local operators with the quantum numbers of spin-1/2 
baryons one can choose from. Here we work with the general form of the nucleon interpolating field parameterized as 
follows: 



m = 2e 



abc \ ^(„,aT ( 



t=l 



K^ {x)CJ[d\x))4u\x), 



(3) 



with JI = I , Ji = J2 ~ 75 and J2 — P, which is an arbitrary parameter that fixes the mixing of two local operators. 
We would like to note that when the choice /3 = — 1 is made the interpolating field above gives what is known as loffe 
current for nucleon. Here u{x) and d{x) denote the u- and d-quark fields, respectively, a, 6, c are the color indices 
and C denotes charge conjugation. 

The short-distance physics corresponding to high momenta p'^ and q^ is calculated in terms of quark and giuon 
degrees of freedom. Inserting the interpolating field in Eq. ^ into the correlation function in Eq. ^ , we obtain 



1 /■ ^ 



+ (J|)p„ [{C4fSu{~x)ia 



t^-'i-zp 



Ae- 



abc / 



Kix)4iO)d^O)\N) 
where u and d denote the quark fields. Su,d represent the quark propagator 
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i-K^X^ 



(5) 



Here the first term gives the hard-quark propagator. The second term represents the contributions from the nonper- 
turbative structure of the QCD vacuum, namely, the quark and quark-gluon condensates. These contributions are 
removed by Borel transformations. The last term is due to the correction in the background gluon field and gives 
rise to four-particle (N = 4 with an additiona.1 gluon) and five-particle (N = 5 with an additional quark-anti-quark 
pair) baryon distribution amplitudes. In Ref. [22l |. higher Fock states (iV > 3) of the nucleon wavefunction have been 
considered and it has been shown by an application to electromagnetic form factors of the nucleon that the valence 
Fock states dominate while A^ = 4, 5 Fock states give only small contributions. Following the common practice, in 
this work we shall not take into account such contributions, which leaves us with the first term in Eq. ([5]) to consider. 
The matrix elements of the local three-quark operator 



4e 



abc / 



\u'^{aix)u'p{a2x)d'^{a3x)\N) 



(6) 



(oi,2,3 are real numbers in the interval [0, 1] denoting the position of the quarks along the ray pointing from the origin 
to the space-time point x) can be expanded in terms of DAs using the Lorentz covariance, the spin and the parity of 
the baryon. Based on a conformal expansion using the approximate conformal invariance of the QCD Lagrangian up 
to 1-loop order, the DAs are then decomposed into local nonperturbative parameters, which can be estimated using 
QCD sum rules or fitted so as to reproduce experimental data. 

The long-distance side of the correlation function is obtained using the analyticity of the correlation function, which 
allows us to write the correlation function in terms of a dispersion relation of the form 



<(p,„ = if 



ImH;y(g) 

(.S-p'2) 



ds 



The ground-state hadron contribution is singled out by utilizing the zero-width approximation 

Im n;^ = nS{s m%){0\v^\N{p')){N{p')\T,,\N{p)) 

expressing the correlation function as a sharp resonance plus continuum with the continuum threshold, sq. The matrix 
element of the interpolating current between the vacuum and baryon state is defined as 

{0\r,''\N{p,s))^XNvip,s) 

where Xn is the nucleon overlap amplitude and v{p, s) is the nucleon spinor. 

The QCD sum rules are obtained by matching the short-distance calculation of the correlation function with the 
long-distance calculation. Using the most general decomposition of the matrix element in Eq. (jB]) we obtain the 
following sum rules for the tensor form factors: 

Ht =4^e™"/^^^ |m^y(l, 1) + ^ ^y (2, 4* + 1) + ^y (1, 13) + J2 ^nZ{^i + 1) + m%xoZ{13)\ , 

Et =^e™«/^'^^ I m^y(l, 1) + ^ ^y (2, 4* + 1) + ^y (1, 13) + J2 ^^Z^^^ + 1) + '«^^o^(13)| , (7) 

1 f 2 3 32 ^ 

ilT =^e'"-/^^' <^ mwy(l, 1) -K ^ 9^^(2, 4* + 1) + ^1^(1' l^) + E '^^^(^z + 1) + m^xoZ(13) \ , 

I, Z— 1 2 — 1 ) 

where we define 

y(fc, n)= \ -^e-^(*^)/*^' {[(1 - mn + (1 + /3)^„+i] + [(1 - /3)F„+2 + (1 + /3)^„+3]} , 

^H = n2^ 2 2 [(1 - mFn + ^„+2) + (1 + P){Fn+l + F^+s)] , 



(8) 



for Ht form factor at the structure i^a^j^^, for Et form factor at the structure 7^9;^ — 7^9/^ ^^nd for iJy form factor at 
the structure Pf^qi, — Puq^i- The explicit form of the functions that appear in the above sum rules are given in terms 
of DAs as follows: 
for Ht'- 



Fi=2 f ' dh i-Pi +T1-T2 + Tj] ih,t2, 1-h- ts), 
Jo 

F2= dh h^i +A2+ 2A3 + V1-V2+ 2^3] (ti, t2, 1 - tl - t2), 

Fz = 2 dh [-Ai +A3-V1- V3] ih, 1 - <i - t2,t2), 

Jo 

F4 = 2 / dh i-Pi + Si- 2Ti -T3- Tj] {h,l-h- t2M). 
Jo 



10 
F5 = 0, 
Fe = 0, 

f-l-t2 



Fr= f ' dh [^f + 1^/'] (ti, 1 -h- t2M), 
Jo 

(•l-t2 

Fg = / dh [iT^'] (tl, 1 - tl - t2,t2), 

Jo 

Fg = -4 [' dX [ dp [ dh [Ti -T2-T5 + n- 2Tr ~ 2Ts] {h,P, l-h- p), 
Jl Ji Jo 

/■*2 r^ /■! — p 

Fio = 4 / dX dp dh [Ai -A2 + A3 + A4-A5 + Ae-Vi+V2 + V3 + V4 + V5~ Ve] {ti,p, 1 - ti - p), 

Jl Jl Jo 

pt2 pX f^ — p 

Fn=/ dX dp dh[Ai-A2 + A3+A4-A5+Ae + Vi-V2-V3-Vi-V5 + V6]{h,l-h-p,p), 

Jl Jl Jo 



Fi2 = - f dX f dp f dh i-Ti + 3T2 - 2T3 - 2T4 + 3T5 -n + ATj + ATs] {ti,l-h- p, p), 
Ji Ji Jo 

Fi3 - 0, 
Fu = 0, 
Fi5 = 0, 
^16 = 0, 



for Et'- 



Fi=2 dh [-Pi + Si + 2Ti -T3- Tj] (fi, t2, l-h- ia), 



F2 = 2 



dh[-Ai + Vi] {ti,t2,l - h - t2), 



F3 = 2 / dh [-Ai +A2-V1+ V2] [h: 1 - <1 - t2,t2), 



Fa 



4/ dh[Tl]{h,l~ti-t2,t2), 

JO 



l-t2 



F5 = -2 / dii [T^'] ih,t2, l-h- t2), 

Jo 

rl-t2 

Fe^-2 dh [-Af + V,^'] {h,t2, l~h- t2) 



Ft^O, 



Fs=4.f 'd<i[Tf](ti,l-ti-t2,t2), 
"'0 

Fg = 2 / d\ [ dp I dh [-Ti +T2+T5-Te + 2Tr + 2X5] {h,P, l-h- p), 
Ji Ji Jo 

Fio = 0, 

Fii=2 [' d\ [ dp j dh [Ai - A2 + A3 + Ai - A^ + A^ + Vi ~V2 -V3 -Vi -V^ + V^] {hA-h- P, p), 
Ji Jl Jo 

F12 = -A [' d\ [ dp j dh [T2 -T3-Ti + T5 + Tj + Ts] {h, l-h- P,p), 
Jl Jl Jo 

Fi3 = 2 / dh / dp[Pi-P2 + Si-S2-2Ti-2T2 + T3+TA + 2n + 3TT-Ts]{h,l-h-p,p), 

Jt2 Jo 

j-1 fl-tl 

F14--2 / dh / dp[Ai+A2 + 2A3-Ai-2A,, + Vi + V2-2V3 + VA-2V5]ihA-h-p,p), 

Jt2 Jo 

Fi5 = -2/ dp dh[Ai-A2 + A3-Vi + V2 + V3]{hA-h-p,p), 

Jt2 Jo 

F16 = 2 / dp j dh [Pi -P2-S1 + S2 + 2Ti -T3-Ti-Tj- Ts] ih,l-h- P, p), 

Jt2 Jo 



and for Ht'- 
Fi^O 
F, 
F3=0 



l-t2 

2= / dh[Vl-Ai]{hMA-tl-t2) 





rl-t2 



/•L-Z2 

Fi^-2 dti[Ti]{ti,l-ti-t2,t2), 

Jo 

^5=0, 

F6 = - f ' dh [Vl' - ^f ] (tl,i2, l-h- t2), 

Jo 

Fs=2 dii[T/^](ii,l-ii-i2,i2), 

Jo 



F9=0, 
Fio = 0, 
Fii = 0, 

Ct2 p\ pl-p 



Fi2 = -2 [' dX [ dp j dh [T2 -T3-T4 + T5 + T7 + Tg] (ii, 1 - fi - p, p), 
Ji Ji Jo 

Fi3 = / dh dp [Pi -P2 + S1-S2- 4Ti - 2T2 + 3T3 + Ti + 2T5 + 5T7 - Tg] (h, l-h- p,p), 

Jt2 Jo 

Fu = / dti dp [Ai +A2 + 2A3 - A4 - 2A5 + V1 + V2- 2^3 + 14 - 2^5] {ti,p, l-h- p), 

Jt2 Jo 

/■I jl-P 

Fi5 = 2 dp dti[Ai-A2 + A3-Vi + V2 + V3]{tul-ti-p,p), 

Jt2 Jo 

Fie = - f dp f dh [Pi -P2-S1 + S2 + 2Ti -n-Ti-Tr- Tg] [hA-h- p, p). 

Jt2 Jo 

The DA's Pi, Ai, Si and Ti appear in the most general decomposition of the nucleon matrix element [23| and the twist-5 
DA's Ai^ , y/^and Tf^ have been calculated in Ref. [20|. In obtaining the sum rules, Borel transformation is applied to 
eliminate the subtraction terms in the spectral representation of the correlation function and to exponentially suppress 
the contributions from excited and continuum states. The contributions of the higher states and the continuum are 
modeled using the quark-hadron duality and subtracted. Both of the Borel transformation and the subtraction of the 
higher states are carried out using the following substitution rules: 

d.-J^^^-C^-^p{x)e-<^y'^\ 
[q-xpY J^^ X 



[q-xpY ' A'P J.^x^'^''^'" ' Q-^+xlmj 



I xo -^ ^ ^ -^0 N 

where 

\ — X 
s[x) = (1 — a;)m« H Q^. 

X 

M is the Borel mass and xq is the solution of the quadratic equation for s = sq: 



Xo 



^(Q2 + So - m%)^ + im%{Q^) - (Q^ + sq - m%) 



/(2"4), 



where So is the continuum threshold. 

To obtain a numerical prediction for the form factors, the residue, Atv, is also required. The residue can be obtained 
from the mass sum rule, which is given by: 

2 



where x = sa/M , and 
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FIG. 1. Tensor form factors Ht{Q^), Et{Q^) and Ht{Q^) as a function of Q^ above Q^ = 1 GeV^. We also show our fitted 
form factors (see the text for details) and give the results for Ht as obtained from xQSM [1Q| and from lattice QCD [I3 



0.8 GcV^ and itin ^ 0.94 GeV. 



Wc use the following parameter values: {qq) = —(0.243)'^ GeV"^, n 

The DAs of the nucleon are given in Ref. [20[ and they are expressed in terms of some nonperturbative parameters 
which are calculated using QCD sum rules or phenomenological models. In this work, we use Chernyak-Zhitnitsky-like 
model of the DAs (see Ref. [20] for details) to calculate the tensor form factors. Note that unlike other form factors 
such as axial ones, the tensor form factors are renormalization-scale, /i, dependent llj . Our results correspond to 
fj^ = 1 GeV^ as we are using the numerical values of the DAs at this scale [2^ . 

There are several parameters in the sum rules that need to be determined. We use the continuum threshold value 
So ~ 2.25 GeV^, which is pretty much fixed in the literature from nucleon spectrum analysis. We also give our 
numerical results at Af^ =2 GeV^. The sum rules and hadron properties are expected to be independent of the 
mixing parameter /?. In order to see if we can achieve such an independence, we make a reparameterization using 
j3 = tan 9 and observe that a stability region with respect to a change in the mixing parameter can be found around 
cos 9 ~0, which we use to produce our results. 

We plot the tensor form factors Ht{Q^), Et{Q^) and Ht{Q^) as a function of Q^ in Fig. [T] We show the form 
factors above Q^ = 1 GeV^. We also give the results for Ht as obtained from xQSM jig and from lattice QCD 17| . 
which have been calculated at a renormalization scale of ji^ = 0.36 GeV^ and ji^ — A GeV^, respectively. We use the 
following relation [2^ to relate these results of the form factors to those at /i^ = 1 GeV'^: 



x{^,^) 



asip-"^ 



4/27 r 



where fii is the initial renormalization scale and 



1- 



337 

4687r 



[asinf) -Usin'^)] 



X{p!), 



An 



91n(MVA|cD) 



64lnln(^VA|cp) 
81 ln(/iVA 



QCD^ 



(10) 



(11) 



with Aqcd = 0.250 GeV and N^ = Nf = 3. 

In order to extrapolate the form factors to low-momentum region, we have fitted QCDSR results to various forms 
including monopole and dipole, which unfortunately fail to give a reasonable description of data with a two-parameter 
fit. On the other hand, an exponential form, viz., 



(12) 

larger than a 
2 GeV^. Our 



^t(Q') = FT(0)exp[-gVm|,], 

allows a plausible description of data with a two-parameter fit. Note that the LCSR are reliable for Q^ 
few GeV^. We have, however, tried three fit regions, namely, Q^ > I GeV^, Q^ > 1.5 GcV^ and Q^ > 
results are given in Table HI EFtJ O) values are to be compared with those from experiment as Ht{0) - 
from xQSM as HriO) = 1.23 [3, from external-field QCD sum rules method as Ht{0) = 1.29 ± 0.51 [12|, and from 
lattice QCD as Ht{0) = 1-15 ± 0.15 [131 at /x^ = 1 GeV^. We find that our results for Ht agree well with those from 
other theoretical approaches, when we make a fit from a more reliable high-Q^ region, namely from Q^ > 2 GeV^ (all 
theoretical methods overestimate Ht when compared to the phenomenological result in [3]). In Fig. [1] we also show 
our fitted form factor as well as those from other theoretical approaches. In the high-Q^ region we observe that our 
Ht results compare well with lattice-QCD results. With the finding that an exponential fit above 2 GeV^ gives more 



TABLE I. The values of exponential fit parameters, namely Ft{0) and rriT, of tensor form factors. We give the results of the 
fits from three regions. 



Form Factor Fit Region (GeV^) Ft(0) mr (GeV) 



Ht 



Et 



Ht 



[2.0-10] 


1.15 


1.35 


[1.5-10] 


1.52 


1.25 


[1.0-10] 


2.11 


1.13 


[2.0-10] 


0.96 


1.11 


[1.5-10] 


1.33 


1.03 


[1.0-10] 


1.92 


0.94 


[2.0-10] 


0.43 


1.10 


[1.5-10] 


0.63 


1.01 


[1.0-10] 


0.97 


0.91 



accurate results for Ht, we quote our final results for the other two form factors at Q^ = as extrapolated from this 
region as Et{0) = 0.96 and Ht{0) = 0.43. 

In conclusion, we have extracted the isovector tensor form factors of nucleon, which play an important role in 
understanding the transverse spin structure of the nucleon when related to the quark helicity-flip generalized parton 
distributions via their first moments. To this end, we have employed LCSR method to leading order in QCD and 
included distribution amplitudes up to twist 6. We have found that our results for the Ht form factor compare well 
with lattice-QCD results in the high-Q^ region and with those from other theoretical approaches when extrapolated 
to low-Q^ regions via an exponential fit. We have also calculated the tensor form factors Et and Ht for the first time 
in the literature. We will extend our work with a detailed analysis of tensor form factors with respect to each quark 
contribution so as to extract isoscalar, octet and singlet form factors as well as those of strange baryons, in a future 
publication. 
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